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ABSTRACT 



Aims. Our aim is to study the photospheric flux distribution of a twisted flux tube that emerges from the solar interior. We also report 
on the eruption of a new flux rope when the emerging tube rises into a pre-existing magnetic field in the corona. 
Methods. To study the evolution, we use 3D numerical simulations by solving the time-dependent and resistive MHD equations. We 
qualitatively compare our numerical results with MDI magnetograms of emerging flux at the solar surface. 

Results. We find that the photospheric magnetic flux distribution consists of two regions of opposite polarities and elongated magnetic 
tails on the two sides of the polarity inversion line (PIL), depending on the azimuthal nature of the emerging field lines and the initial 
field strength of the rising tube. Their shape is progressively deformed due to plasma motions towards the PIL. Our results are in 
qualitative agreement with observational studies of magnetic flux emergence in active regions (ARs). Moreover, if the initial twist 
of the emerging tube is small, the photospheric magnetic field develops an undulating shape and does not possess tails. In all cases, 
we find that a new flux rope is formed above the original axis of the emerging tube that may erupt into the corona, depending on the 
strength of the ambient field. 
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1. Introduction 



\ Active regions are often associated w ith episodes o f magnetic 
flux emergence from the solar interior dZwaanl (Il985l) and refer- 
ences therein). An important question then is, what is the evolu- 
tion of the magnetic field configuration at the photosphere dur- 
ing emergence? Observations of emerging flux regions (EFRs) 
as recorded at the photospheric level, show that they consist of 
. two main flux bundles of opposite magnetic polarity that may be 
the manifestation of an emerging flux tube. There is strong evi- 
dence that, in many EFRs, the rising magnetic fields are twisted. 
The idea of a flux rope configuration has been supported by 
photospheric measurements and observations of emergi ng fields 
in normal and complex (the so-called S sunspot) ARs (|Tanakal 
I1991I: iLites et al.lfl995l: iLeka et al.lfT996l: ICanou et al.ll2009h . 
A common feature in EFRs is the presence of magnetic tongues 
or tails, which are connected with the main p olarities on the two 
sides of the PIL of the AR dLi etalJl2007E ICanou et alJ l2009t 
IChandra et afll2009l) . The appearance of magnetic tails is inter- 
preted as the result of th e emergence of twisted ma gnetic field 
lines at the photosphere (lLopez Fuentes et ai1l2000l) . However, 
a study of how the formation and evolution of the tails depend 
on the physical pro perties of the emerging field is still missing. 
ICanou et al.l (l2009h used SOHO/MDI magnetograms and they 
reported on the existence of tails, which formed along the PIL 
and accompanied the emergence of magnetic flux in the region 
NOAA AR 1808. The shape of the tails was deformed during 
the evolution of the system. They also used the THEMIS vector 
magnetogram to reconstruct the coronal field (via a nonlinear 
force-free model) and found evidence for a pre-emptive twisted 
flux tube above the emerging field. 

In this paper, for the first time, we focus on the photospheric 
distribution of an emerging flux tube and the formation of the 



tails, showing the relationship between the topology of the tails 
and the initial tube parameters. We co mpare some of the nu- 
merical results with the observations bv lCanou et all d2009h and 
we find a preliminary, qualitative agreement. Secondly, we re- 
port on the emergence of the tube into a magnetized corona 
and the subsequent cor onal eruption of a flux rope. Similar 
to previous experiments < 



onal eruption oi a flux rope, bimilai 
dMagaralhoOll: iMancnester et al.ll2004 



[ Archontis & Hood 120081: 1 Archontis & Torokll2008l: iHood et alJ 
120091) . we find that the emerging twisted flux tube and the coro- 
nal rope are two distinct structures. More importantly, we find 
that the evolution of the erupting rope (ejective vs confined erup- 
tion) depends on the strength of the ambient field. 

2. Model 

The results in our experiments are obtained from a 3D MHD 
simulation . The basic setup of th e experiment follows the sim- 
ulation by I Archontis et al.l (120051) and consists of a hydrostatic 
atmosphere and a horizontal twisted magnetic flux tube. All 
variables are made dimensionless by choosing photospheric val- 
ues of density, p p h = 3 x 10~ 7 gcm~ 3 , pressure, p p h = 1.4 x 
10 5 ergs cm -3 , and pressure scale height, = 170 km, and 
by derived units (e.g., magnetic field strength = 1300 G, 
velocity V p h = 6.8 km s" 1 and time f p h = 25 s). The atmo- 
sphere includes a subsurface layer (-10 < z < 20), photosphere 
(20 < z < 30), transition region (30 < z < 40), and corona 
(40 < z < 110). The numerical domain has a dimensionless 
size of [-70,70] x [-80,80] x [-10, 110] in the longitudinal 
(x), transverse (y) and vertical (z) directions, respectively. The 
magnetic flux tube is imposed \AMm below the surface along 
the y-axis. The radius of the tube is 425Km. The axial field is 
defined by B y = Bq exp(-r 2 /7? 2 ) and B<p = arB y , where r is 
the radial distance from the tube axis and a is the twist per unit 
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Fig. 1. Top: SOHO/MDI magnetograms during flux emergence in NOAA AR 10808. Bottom: Magnetograms produced in the 
numerical experiments, at z = 21 and t = 40, t = 125 and t = 175 for the panels Id- If respectively. Arrows show the horizontal 
component of the magnetic field, Bh or (panels Id and le, black) and velocity field, Vh or (panel If, red). 



length. The twist of the fieldlines around the axis of the tube is 
uniform. The tube is made buoyant by applying a density pertur- 
bation Ap = [p t (r) I p(z)] p(z) exp(-y 2 M 2 ), where p t is the pres- 
sure within the flux tube and A is the buoyant part of the tube. 
We perform four experiments: El (B = 3, a = 0.4, A = 10), E2 
(B = 5, a = OA, A = 10), E3 (B = 5,a = OA, A = 20) and 
E4 (B = 5, a - OA, A - 10). In the numerical domain, we use 
periodic boundary conditions in horizontal directions and closed 
boundaries with a wave damping layer in vertical directions. 

3. Results 

3.1. Emergence into the photosphere 

Figure 1 (panels la-lc) shows the evolution of the emerging field 
in the NOAA AR 10808. At the beginning (panel la) there is 
a clear appearance of a bipolar region at the photosphere with 
a North-South orientation. The two polarities progressively di- 
verge from each other in an approximate East- West direction 
(panels lb, lc). During the evolution of the system, two elon- 
gated tails or tongues are formed in the wake of the two polari- 
ties (panel lb). Initially, the tails possess an apparently coherent 
shape but as time goes on their structure appears to be more frag- 
mented on the two sides of the PIL (panel lc). 

Panels Id- If show the photospheric distribution of the 
emerging field in our numerical experiments. Panel Id shows 
the bipolar appearance of the emerging field, shortly after it 
intersects the photosphere. The North-South orientation of the 
bipolar field is due to the strong initial twist of the flux tube. 
Eventually (panel le), the two main polarities drift apart toward 
an East- West orientation. Similar to the observations, they are 
followed by magnetic tails that develop an intricate geometrical 
shape. The projection of the horizontal component of the mag- 
netic field (arrows) is overplotted onto the magnetograms of pan- 
els Id and le. At t = 40, the direction of the horizontal magnetic 
field vectors shows a normal configuration, i.e. from the positive 
to the negative polarity, at the PIL. Later on, as more magnetic 
flux emerges from the solar interior, the direction of the mag- 
netic field reveals a dominant inverse configuration along the 



PIL. This is due to the rise of the original axis of the twisted flux 
tube above this height (z = 21). However, the main axis does not 
emerge above 2-3 pressure scale heights, as has been shown 
in previous experiments of flux emergence (lFanll2001l:lMagaral 
1200 ll:l Archontis et al.ll2005h . 

At a later stage of the evolution (Panel If), the elongated 
tails develop fingers seperated by dips along the curved PIL. In 
the fingers, the magnetic field remains strong (around 70% of the 
maximum value of B z at this height). At the dips, the magnetic 
field is weak and the plasma density is relatively small. In fact, 
we find that there is a good correlation between the location of 
the dips and sites where plasma is moving in the transverse direc- 
tion. There, the converging flows may reach values up to 3Km/s 
and the kinetic energy density becomes larger than the magnetic 
energy of the field. Thus, it seems that the shape of the magnetic 
tails is deformed due to inflows that are able to compress and 
advect the magnetic field. The origin of the inflows depends on 
the evolution of the total pressure (P tot = magnetic + gas pres- 
sure) at photospheric heights. Panel 2a shows the distribution of 
P tot at z = 25, when the outer magnetic field has expanded into 
the corona. Due to the rapid expansion, a total pressure deficit 
has developed at the central area of the EFR and so the plasma 
moves towards the small pressure, and deforms the tails. 

The link between the appearance of the tails and the topology 
of the fieldlines is shown in panel 2b. The yellow fieldlines have 
been traced from a far edge of the fingers (at x - -11, y = 
0). These are the outermost fieldlines with a strongly azimuthal 
nature. The blue fieldlines are highly twisted and are traced from 
the fingers of the tails that are closer to the PIL. They make a 
full turn around the main axis of the emerging tube connecting 
the central area of the two tails. The red fieldlines have been 
traced from the region closer to the main positive polarity of 
the field. They are very weakly twisted, possessing an arch-like 
bundle of fieldlines, joining the two sunspots. These fieldlines do 
not go through the tails. The above configuration shows that the 
appearance of the tails is due to the projection of the azimuthal 
component of the magnetic field at the photosphere. 
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Fig. 2. Top: The colour-scaled maps correspond to the P tot (2a) and B z (2b, 2c). Contours show B z and arrows the horizontal velocity. 
Time is t = 165 and z - 25, all for experiment El. Bottom: Distribution of B z for E2 (2d), E3 (2e) and E4 (2f). 



Panels 2c-2f show the magnetic flux distribution at the pho- 
tosphere for the experiments E1-E4 respectively. We take as a 
reference case the El and we examine the effect of varying the 
initial field strength B, X and a on the appearance of the tails. 
For comparison, we consider the configuration of the field at a 
certain time for all experiments. The increase of B (in E2) re- 
sults in keeping a coherent shape of the tails for a longer time 
period: at t = 165, the tails in E2 are less fragmented than in 
El. This is due to the fact that the total pressure within the EFR 
is large enough for the tails to be distinctively deformed by the 
inflows. However, we should emphasize that the shape of the 
tails is altered at a later time, when the two sunspots have seper- 
ated enough and the magnetic field in the EFR becomes weak. 
The increase of A (panel 2e) affects the downward tension of the 
fieldlines upon the buoyant part of the emerging field. The ten- 
sion is less in the E3 and the field is emerging at the photosphere 
relatively faster. Thus, at a certain time, the magnetic field at the 
photosphere appears stronger in E3 than in E2. As we mentioned 
above, the stronger the magnetic field the less effective is the de- 
formation of the tails' shape. This is clearly shown in panel 2e, 
compared to the E2 (panel 2d). Also, the appearance of the tails 
critically depends on the initial twist of the emerging field. In 
E4, the twist parameter a is equal to 0.1 and the emerging field 
is almost horizontal and parallel to the E-W direction, shortly 
after its arrival to the photosphere. We find that there is no tail 
formation when the emerging field has a < 0.2. In this case, the 
EFR consists of the two sunspots and patches of magnetic flux 
with mixed polarity on the two sides of the PIL. Some of these 
photospheric flux segments are connected with the same field- 
lines, possessing an overall undulating magnetic system. This 
is reminiscent of the "sea-serpent" configuration, which is pro- 
duced during the emergence of a magnetic flux sheet. The latter 



develops undulations when it beco mes unstable to the Parker in- 
stability (lArchontis & Hoodll2009l) . 



3.2. Eruption into the corona 

In Section 3.1, we showed that the photospheric fingerprints of 
the EFR in El consist of features (e.g. tails) with a similar con- 
figuration to observed ARs (e.g. the AR 10808). In addition, the 
activity in the region NO AA AR 10808 is k nown to lead to fila- 
ment and CME eruption (ICanou et al .1120091) . Thus, an important 
question is whether our twisted flux tube model can produce a 
coronal eruption. Our experiment shows that a new flux rope is 
formed above the original axis of the emerging flux tube due to 
reconnection of sheared fieldlines. The reconnection occurs in 
the higher photo sphere/lower chromosphere in a sim ilar manner 
to the model by Ivan Ballegooijen & Marten s (1989). A key is- 
sue is whether this eruption is confined (and, thus, the flux rope 
cannot fully escape into the outer atmosphere) or ejective. In pre- 
vious experiments J Archontis & Toro k ( 2008) found that the in- 
clusion of a pre-existing magnetic field in the corona may induce 
a runaway situation, via reconnection, during which the new flux 
rope fully erupts into the outer solar atmosphere. Here, we per- 
form a similar experiment but using different initial parameters 
for the pre-existing coronal field. Our aim is to study whether 
the field strength of the ambient field affects the rising motion of 
the erupting flux rope. 

The observed magnetogram in the AR 10808 shows that the 
emerging flux is rising into a pre-existing field oriented in the E- 
W direction. The emerging field has a N-S orientation and, thus, 
the relative orientation of the two fields is about 90 degrees. To 
simulate this, we include a horizontal and uniform magnetic field 
in the corona along the y-axis, parallel to the main axis of the 
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Fig. 3. Height-time profiles of the apex of the emerging field 
(solid) and the flux rope (dashed) in experiments B 1 (black), B2 
(red) and B3 (green). 



twisted tube (for example, see I Archontis et al. (2005)). To a first 
approximation, this field may correspond to the upper part of the 
observed AR's field, which is likely to be anchored in the sur- 
rounded diffuse polarities. We find that the field strength of the 
ambient field {B am b) plays a critical role in the eruptive motion 
of the new flux rope. Figure 3 shows the height- time profile of 
the front of the emerging field (solid lines) and the center of the 
new flux rope (dashed lines) for three experiments (Bl, B2 and 
B3) where: B amb = 0.003 (Bl, black lines) , B amb = 0.015 (B2, 
red lines) and B am b = 0.03 (B3, green lines) respectively. The 
heights are calculated at the vertical midplane after the emerg- 
ing field enters the transition region. 

In Bl, the front of the expanding tube rises slowly within 
the magnetized corona and eventually it saturates at a height of 
z ~ 96. The new flux rope is formed at the low atmosphere at 
t ^ 95 and, thereafter, it follows a similar evolution to the enve- 
lope field of the expanding tube. Firstly, it rises almost linearly 
with time but then it reaches an equilibrium where the magnetic 
pressure force is balanced by the tension of the fieldlines. In this 
case, the eruption is confined: the flux rope is trapped within 
the envelope field. In B2, the apex of the emerging field reaches 
lower heights during its rising motion. This is because it comes 
into contact with an ambient field that is stronger and able to 
delay the emergence. At the same time, a considerable amount 
of the rising magnetic flux is removed from the envelope field 
due to reconnection. As a result, the distance between the new 
flux rope and the front of the envelope field is reduced. As more 
magnetic layers above the flux rope are peeled off, the down- 
ward tension of the envelope fieldlines decreases. Eventually, 
the flux rope experiences an ejective eruption reaching the up- 
per boundary of the domain very quickly. Due to the short dis- 
tance between the erupting rope and the closed top boundary, the 
velocity of the center of the flux rope is restricted to 197 Km/ s. 
However, the plasma underneath the flux rope is rising with even 
higher velocity at « 350Km/ sec. This is a reconnection jet that 
is formed due to internal (i.e. within the EFR) reconnection of 
fieldlines and helps the flux rope to accelerate during its erup- 
tion. According to these calculations, it is possible that the rise 
of the flux rope might account for a CME-like eruption. 



In B3, the eruption of the flux rope is triggered earlier. Again, 
this is because the stronger ambient field reconnects more effec- 
tively with the flux above the rope and removes more magnetic 
layers from the emerging system. However, for the same rea- 
son, the front of the envelope field rises with a slower rate and 
the distance between the new flux rope and the front decreases. 
As a result, soon after the triggering of the ejective eruption, the 
erupting rope collides with the front and loses its distinct circular 
shape, possibly due to reconnection with the ambient field. After 
the collision, the leading edge of the emerging system is lifted 
up for a few pressure scale heights. However, it does not recon- 
nect effectively with the magnetic flux above it, and eventually 
reaches a quasi-static state at a height of z ~ 80. Thus, in B3, the 
ejective flux rope is trapped by the dominant ambient field and 
not by the envelope field. 

4. Summary and Discussion 

In this paper, we have presented a 3D model to study the emer- 
gence of a twisted flux tube throughout the solar atmosphere. 
Our model gives new insights into the photospheric distribution 
of the emerging magnetic field: it consists of a bipolar region 
and tails on the two sides of the PIL. The appearance of tails re- 
veal that the emerging magnetic field is twisted. For small twist, 
the emerging field possess undulations. Our results predict that 
the irregular structure of the tails is due to the interplay between 
the flows and the dynamical evolution of the magnetic field. The 
configuration of the emerging field a t the photosphere i s in qual- 
itative agreement with observations (iCanou et al.ll2009l) . 

In agreement with previous simulations, our experiments 
show the eruption of a flux rope, which is formed above the orig- 
inal axis of the emerging tube. For the first time, we find that the 
field strength of a pre-existing coronal magnetic field is a crucial 
parameter affecting the eruptive phase of the rope. Under the 
specific conditions of the present experiments, we found that the 
eruption is ejective when 0.01 < B am h < 0.02. For other values, 
the eruption is confined within the envelope field. 

The aim of these experiments is not a direct comparison with 
the observations, but rather to suggest possible mechanisms that 
drive the dynamical behaviour of the system. Further experi- 
ments are required to verify the effect of the initial parameters 
(e.g. field strength, radius, initial atmospheric height and twist, 
etc.) of the twisted flux tube and the pre-existing field on (a) 
the characteristics of its photospheric appearance (formation and 
evolution of the tails, shear and transverse flows, etc.) and (b) the 
dynamics of the associated eruption. 
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